x=[0 0.005 0.01 0.05 0.1 0.2 0.5 0.7 0.9 0.95 0.99 0.995 1];
t=[0 0.005 0.01 0.05 0.1 0.5 1 1.5 2];

x=0:0.05:1;
t=0:0.05:2;
m=0;

sol=pdepe(m,@f1,@f2,@f3,x,t);
u1=sol(:,:,1);
u2=sol(:,:,2);

figure;
surf(x,t,u1)
title('u1(x,t)')
xlabel('Distance x')
ylabel('Time t')

figure;
surf(x,t,u2)
title('u2(x,t)')
xlabel('Distance x')
ylabel('Time t')
